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The light scattering problem for a confocal multilayered spheroid has 
been solved by the extended boundary condition method (EBCM) with a 
corresponding spheroidal basis. The solution preserves the advantages of the 
approach applied previously to homogeneous and core-mantle spheroids, i.e. 
the separation of the radiation fields into two parts and a special choice of 
scalar potentials for each of the parts. The method is known to be useful in a 
wide range of the particle parameters. It is particularly efficient for strongly 
prolate and oblate spheroids. Numerical tests are described. Illustrative 
calculations have shown that the extinction factors to converge to average 
values with a growing number of layers and how the extinction vary with a 
growth of particle porosity, (c) 2011 Optical Society of America 

OCIS codes: 290.2200, 290.5825, 290.5850 
1. Introduction 

A detailed knowledge of the optics of inhomogeneous (layered) non-spherical particles is re- 
quired in many scientific and industrial applications. Numerical treatment of these particles 
is a very complicated problem, especially when the particle size is not small and (or) the 
particle shape appreciably deviates from spherical (see [1-3] for a review of available meth- 
ods). However, for one of the simplest cases, multilayered spheroids, rather fast calculations 
of a high accuracy can be performed by the method of separation of variables (SVM) and 
the extended boundary condition method (EBCM). Both methods can be used with spheri- 
cal or spheroidal basis so that the electromagnetic fields are expanded in terms of spherical 
or spheroidal wave functions, respectively [3,4]. As a result, characteristics of scattered ra- 
diation can be calculated by using the same expressions. Note that the methods differ in 
formulation of the boundary conditions (see [4] for detailed discussion). However, the use of 
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the spherical basis is not appropriate for particles of large eccentricity (with aspect ratios 
a/6 ^ 1.5 — 2) that is why one needs to apply a spheroidal basis for these particles when 
geometry of the problem is sorely taken into account. 

The first attempt to develop a solution for multilayered confocal spheroids by SVM with a 
spheroidal basis was made in [5] by using a recursive procedure when passing from one layer 
to the next. The paper did not contain calculations because they require solving a complex 
nonlinear matrix equation for the unknown expansion coefficients. Later the algorithm was 
modified by using the ideas presented in [6] and some numerical results were published in [7] 
for small particles with large refractive indices. Using the SVM approach, an exact solution 
for spheroids with non-confocal layers was obtained but the calculations were published for 
core-mantle particles only [8]. 

Layered axisymmetric particles (including spheroids) were also treated by the T-matrix 
method (e.g., [9,10]) and generalized multipole technique or null-field method (e.g., [11]). 
However, these methods did not provide the adequate numerical results for strongly non- 
spherical particles with a large number of layers (see discussion in [3]). 

In this paper, we consider the scattering of an arbitrary polarized plane wave by confocal 
multilayered spheroids. We have developed the recursive EBCM solution with a spheroidal 
basis suggested in [6] by taking into account inaccuracies found during our numerical real- 
ization of the algorithm (see also [12]). Our solution is based on a special choice of scalar 
potentials which for any next layer can be found by using the potentials of the previous 
layer, the procedure starting from the particle core. These potentials are expanded in terms 
of spheroidal wave functions. The unknown expansion coefficients of scattered radiation po- 
tentials are determined by solving the systems of linear matrix equations. It is important 
to emphasize that the dimension of these systems for layered spheroids does not increase as 
compared to that for homogeneous spheroids which is contrary to the SVM (see, e.g., [13]). 
Calculations show that the method suggested in this paper gives results of high accuracy 
and can be used in a wide range of the particle parameters. This confirms the conclusion 
made in [3] that the EBCM with a corresponding spheroidal basis is more preferable in the 
treatment of multilayered confocal spheroids. 

2. Formulation of the Problem 

The problem of electromagnetic light scattering by a multilayered spheroidal particle is solved 
in the prolate and oblate spheroidal coordinate systems (^, t], (p) which are connected with 
the Cartesian system (x, y, z) in the following way [14,15]: 
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where / = 1, ^ G [l,C)o), rj G [-1,1], G [0,27r) for prolate coordinates and / = —1, ^ G 
[0, oo), ?7 G [-1, 1], G [0, 27r) for oblate coordinates; d is the focal distance. We assume that 
the particle is confocal. This means that the surfaces of layers coincide with the coordinate 
surfaces, and their equations can be written as 

e = 0, (2) 

where j = 1,2,...,N [N > 3 is the number of layers, j = 1 for the outermost surface, 
i.e. the particle boundary, and j = N for the boundary of the core). For such particles, the 
major and minor semiaxes of the shell spheroids, aj and bj, satisfy the following conditions: 

ai-bl = al-bl= ... =a%-b% = (0'. (3) 

Let the time- dependent part of the electromagnetic field be exp {—iut) and E, H be the 
vectors of the electric and magnetic fields, respectively. The vectors E^^\ H^^'^ correspond 
to the field of the incident radiation, E^^\ H^^'^ to the field of the scattered radiation, E^'^\ 
i?*^^) to the field inside the outermost layer, . . ., E^^\ H^^^ to the field inside the (j — l)th 
layer, . . ., E^^^'^\ to the field the particle core. 

We consider a plane electromagnetic wave with an arbitrary polarization propagating at 
an incident angle a to the rotational axis of the spheroid (or the z-axis; see Fig. 1). This 
wave can be represented as a superposition of two components (the magnetic fields can be 
obtained from the electrical ones by using Maxwell's equations): 

(a) TE mode 

E^^'^ = —iy exp [iki {x sin a + z cos a)] ; (4) 

(b) TM mode 

E^^'' = {ix cos a — iz sin a) exp [iki {x sin a + z cos a)] . (5) 

Here i^, iy, iz are the unit vectors in the Cartesian coordinate system, ki = y/eiJTiko is the 
wave number in a medium with the complex permettivity £j and the magnetic permeability 
/ij, ko = 27r/Ao and ki = 27r/Ai are the wave numbers in vacuum and the medium outside 
the particle, respectively. 

As it has been previously shown in [16] (see also [13, 17]), in the case of axisymmetric 
particles, the scattering problem can be solved independently for each term of the Fourier 
expansion of the vectors E^^^ and H^^^ in terms of the azimuthal angle if. In the following, 
we represent all electromagnetic fields as 

E«=eP+4*), #» = i?f)+i?«, ^ = 0,l,2,...,iV+l (6) 



so that E^^ and are independent of the azimuthal angle (p (the zeroth term of the Fourier 
series), whereas the averaging of E^'^ and #f over (f gives zero. Below, the axisymmetric 
problem for the fields e[^\ Hf' and the non-axisymmetric problem for the fields -^2*^ -^2*'' 
are solved independently of one another. 



3. Solution to the Axisymmetric Problem 

Let us consider the scalar potentials 

p« = e;2cos(/^, q« 



Hl^ cos if, 



(7) 



where if}^ are the (/^-components of the vectors ^| , (z = 0, 1, 2, . . . , + 1). If we 
remove the cos(f factor, these potentials coincide with the Abraham potentials within the 
factor h„ 



— {d/2)y (,^2 _ _ j^2^ ^iQj_ follows from Maxwell's equations that the scalar 
potentials satisfy the wave (Helmholtz) equations 



AP(*) + P» = 0, AQ(*) + Q 



ii) 



0. 



The remaining components of the electromagnetic fields e[^^, H^^ can be expressed in terms 
of their azimuthal components. Note that the axisymmetric problem is solved independently 
for potentials V and Q, i.e., for the TE and TM waves. 

In the case of the TE mode (see Eq. (4)), the boundary conditions (the continuity of the 
tangential components of the electromagnetic fields at the interfaces) should be rewritten as 



p(0) +p(l) = -p(2)^ 



p{j) ^ -pO+i)^ 



d 



^2 _ J-p(2) 



(9) 



«=6 



d 



+1) 



(10) 



where j = 2, 3, . . . , A. 

Let us next formulate the problem in the form of surface integral equations. We can 
represent the potential P'--'^ of the radiation in the (j — l)th shell of the particle as (j = 



2,3, 



A + 11 



-pij) = -pU) + p 



>(i) 

B ' 



(11) 

where V^"* has no singularities in the region -Dj_i (and hence in the region Dj enclosed by 
the surface Sj) and the potential "Pg ■* satisfies the radiation condition at infinity. Note that 
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inside the core = v[^^^\ that is, V^^^^^ = 0. Within the framework of EBCM, we 



obtain the system of surface integral equations (see [6] for more details) 



^1 _p{2)(-) 



(,,.-) ^ 

G\ \ d?7'd<y?' = 




>/o 



(12) 



where ^' = 



2(e/-/) 



27r /•TT 



^0 



p(i+i) 



I2j aP(^ + l)(f') 



- f) 



Gj > dr] d(p 
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where ^' = ^ 



exp ikj\f— f'\ 



AiT\f— f'\ 

is the Green function of the wave equation with the wave number kj, j = 2,3, . . . , N . 
The scalar potentials can be expanded in terms of the spheroidal functions [15] 



(13) 
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(16) 
(17) 
(18) 



where j = 2, 3, . . . , + 1. For the incident radiation we obtain the following coefficients 
[13,17]: 

(a) TE mode (see Eq. (4)) 



,(0) 



-2i'iV-2(ci)5i,(ci,cosa), b\''> = 0; 



(0) 



(b) TM mode (see Eq. (5)) 



,(0) 



0, bP = 2i\l^N-,\cr)Su{c^,cosa). 

V /^i 



(19) 



(20) 
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Here, R^mi^^\'^j^O the prolate radial spheroidal functions of the first and third kinds, 
Smi{C'ji v) the prolate angular functions with the normalization coefficients Nmi{cj) [15], and 
the parameter Cj = kj{d/2). 

For the expansion of the Green function in terms of the spheroidal functions, we have [15]: 



oo oo 



G{kj,r,r') 



271 



m=0 l=m 

>^Smi{cj,r]) Smi{cj,ri') cosm(v? - (p'), 



(21) 



■>0m 



where 

1, m = 0, 
0, m ^ 0, 

and ^< = min(^,^')> ^> = max(^,^')- 

We substitute Eqs. (15)-(18), and (21) into the integral equations (12), (13). Taking into 
account orthogonality of the angular spheroidal functions Smi{cj, rf) cos nup on the surface of 
any spheroid, we obtain the linear algebraic equations for the unknown expansion coefficients 
of the potentials considered. In the matrix notation, they have the following form (j = 
1,2,. ..,iV): 



•^31 -^33 



11 



13 



z 



z 



(22) 



B 



where 



A^ 



(i) 

31 



7^[3l(c„e,) A«(a 



/ij+i 



A«(c„c,+i) 7^W(c,+l,e,) 

^^''(9+i,ej,e,+i), 



(23) 



33 



/^i+i 



- 1 



~ A«(c„c,+i) 



^'"(c,+i,0,e,+i), 



(24) 



A 



(j) 
11 



7^W(c„e,) AW(c„c,+i) - ^AW(c„c,+i) 7^W(c,+l,e,) 



- 1 



~ A«(c„c,+i) 



pW(c,+i,e„0+i), 



(25) 



A 



(j) 

13 



7^W(c„e,) A«(c„c,+i) - ^A(^)(c,,c,+i) 7^Pl(c,+l,e,) 

/ij+i 



/ij+i 







~ A«(c„c,+i) 



p[=^i(c,+i,o,e,+i)- 



(26) 
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Above, we introduce the vectors specified by 



°)i?S;)(ci,ei)iViKci) 

:^)/2«(c,,e,)iVi,(c,) 

^■)i?Sf(c„0)iVi/(c,) 



1 

oo 



1 

oo 



where j = 2, 



N and the diagonal matrices 



>v, = -[7^[3l(c„o)-7^W(9,0)]■ 



^ ^Cji ij-ii ij 



(i), 



nl 



(27) 



(28) 



(29) 

(30) 
(31) 



The matrix elements /s!^'^\cj, Cj+i) = < 5^^\cj, Cj+i) [ are integrals of the products of the 
angular spheroidal functions [13]. To derive Eq. (30), we use expression for the Wronskian 



of the radial spheroidal functions [15]. Since z- 



-iN+l) 



B 



0, the system of equations (22) can be 



easily solved relative to the expansion coefficients of the scattered radiation potential 



A2A\ 



(-1) 



A ' 



(32) 



where the coefficients of the incident radiation are given by Eq. (19). The matrices Ai and 
A2 satisfy the relation 



A, 
A2 



-A 
■^11 



(1) 

31 



-^33 



•^31 -^33 



■^11 



^0) 



13 



-A 



(N~l) 
31 



•^11 



(7V-1) 



-^33 
■^13 



•^31 

■^11 



■ (33) 



An important point is that the representation of Eqs. (32), (33) in the recursive form 
needs only one matrix inversion in contrast with the T-matrix representation that requires 
inversions for each layer and one more at the last step (see discussion in [1]). 

For TM mode, the transformation of the above equations for potentials Q^^^ is performed 
by the replacements — ?■ e^, Sj — >■ yU.j, \ and cp"* d!f \ In order to obtain the 

corresponding systems for oblate spheroid one must use the standard replacements c — )■ — ic, 
^ — 2^ and oblate spheroidal functions instead of the prolate ones. For example, in the case 
of oblate spheroids and TM mode, Eq. (23) can be written as: 



A 



(i) 

31 



V}^\-iCj,iij) A(^)(-zCj, -icj+i) 



Sj+l 



- 1 



^i+i 



m - f 



. A«(- 



-tc 



-ICj+i) 
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4. Solution to the Non-Axisymmetric Problem 



The second terms in Eqs. (6) can be represented in the following form: 
(a) TE mode 



E. 



H. 



V X V X ( + V^'^r 



(34) 



(b) TM mode 



p(0 



E. 



■^—V X V X (u^% + 1/»f) , 



(i) 



V X (u^% + , 



(35) 



where the scalar potentials t/^*-* and V^^^ satisfy the Helmholtz equations (8). 

In the case of TE mode (see Eq. (4)), the boundary conditions for scalar potentials have 
the form 



d 



d 
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\ — rf d 



cj J ^2 _ fdr] 



(36) 



where j = 2, 3, . . . , A^. The boundary conditions for the potentials U^'^\ U^^^ and V^'^\ V^^^ 
can be written in a similar way to the axisymmetric part (see Eqs. (9), (10)). 

As in the case of the axisymmetric part, we can derive the integral equations for the scalar 
potentials U'-^^ and V^^^ 



ai - f) 
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/2 



J' 



f e R^\Dj, 



(37) 
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■2 d 



e|-/V2 



+ 



,/2 



e| - N 

reDj 

fe R^\Dj, 



(38) 
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terms of the 



f/(0) 


oo oo 
m=l /=m 


(0) 
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d{1) 
^m« 


(ci,0 Sr,uici,r]) 


COS rmp, 


(39) 


[/(I) 
1^(1) 


CXI oo 

= EE 

m=l Z=m 


"mi 

"m; 


p(3) 
^mi 


(ci,0 'S'w(ci,?7) 


cosmif, 


(40) 




oo oo 
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"m« 
"ml 


d(1) 
Rml 




cos rmp, 


(41) 


TjU) 
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oo oo 
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■-mi 




cos rrnp, 


(42) 



where j = 2, 3, . . . , iV + 1. For th 

to (see [13,17]) 



WliClC J — ^, . . . 

radiation are equal 

For TM mode, the coefficients have the opposite sign and the muhiphcand ^J Sij ^\ 
I. r35lV 



Eqs. (5), (35)). 



(43) 
(see 
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Substituting Eqs. (39)-(42), and (21) into integral equations (37) and (38), we obtain 
infinite systems relative to the unknown expansion coefficients. The systems can be written 
in the matrix form [j = 1,2,..., A^) 



^1 



z 



where the vectors and matrices have the block structure 



^A 
^B 



kic'i]Rt]{c,,QN^i{c,) 



A. 



U) 

ik 



oi(j) ai(i) 
^ik,B ^ik,B' 



(44) 



(45) 



(46) 



21 



(i) 

31,A 



{7^[^kc„e.) A(-)(c„c,+i) - ^A(-)(c„c,+i) 7^W(c,+l,0) 



m). 



+ 



+ 



-g(™) (c„ c,+i, 0) (/ - 2e|g('")(c,+i, c,+i, Q)] 



;pW(c,+i,0,0+i), 



^Q(™)(c„9+i,0)i?^"nci+i,c,+i) 



(47) 



(i) 

31,A 



+ 



Jlj_ 



1 /e,Q^"Hc„c,+l,e,)^("^)(c,+l,c,+07^™(9+l'e,) 



X 



(48) 
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(j) 

31,B 



e 



JXj_ 



(49) 



(i) 

31, B 



+ 



g('")(c„c,+i,o)(/ - 2e|g('")(c,+i,c,+i,o))] 



(50) 
(51) 

Q^'^Hc,-, 9+1, 0) = A("^)(9, 9+i) Q(™)(9+i, 9+1, 0) (52) 

and / = {5nz}m ^'^^ ^'^^^ matrix. Here, m is the azimuthal index that runs from unity to 
infinity. The subscripts and superscripts of the matrices have the same meaning as in the 
previous section. The elements of the remaining matrices in Eq. (46) can be obtained from 
Eqs. (47)-(50) as explained above (see Eqs. (23)-(26)). The matrix elements r'^™)(9,9) = 
{ln:\cj,c,)r, K^'^Kc.^c,) = {47^(9, 9)}~ and i?(™)(9,9) = {^17^(9, are the 

integrals of products of the angular spheroidal functions and their derivatives (see [13, 17]). 

The system (44) can be easily solved for the expansion coefficients of the scattered radiation 
potential (cf. Eq. (22)) 



7(1) 



(-1) 7(1) 



where 



?(i) 



^i«S^Sl(ci,6)iVw(ci) 




m 
00 



(53) 



(54) 



(55) 
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and the matrices Ai and A2 satisfy Eq.(33) but have the block structure (see Eq. (46)). 

For the TM mode, we can obtain infinite systems for the unknown expansion coefficients 
of the scalar potentials by replacing fij — )■ ej, Sj — fij in the above relations. Note that in 
the case of fij = 1, their form is much simpler than in the corresponding case of the TE 
mode (see also [13,17]). 

5. Characteristics of Scattered Radiation 

Using the expansion coefficients of the scattered field for the TE and TM polarizations, we 
can calculate elements of the scattering matrix (see the corresponding expressions in [13]) 
and the integral characteristics of the scattered radiation (e.g., the cross sections of a particle 
for extinction Cext, scattering Csca, absorption Cabs, and radiation pressure Cpr). These cross 
sections are products of the corresponding efficiency factors Q and the viewing geometric 
cross section of a spheroid (the area of the particle shadow) 

C = GQ, 

where 

1 /2 

G{a) = nbi (a^ sin^ a + bl cos^ a) for prolate spheroids, (56) 

1 /2 

G{a) = nai (a^ cos^ a + sin^ a) for oblate spheroids (57) 

and Oi and bi are the major and minor semiaxes of a multilayered spheroid. The efficiency 
factors for extinction can be found as 

i^'' ai^^ Sii{ci, cos a) 



4 

Qext ~ _, -|^^2 



1=1 



c?[(e?-/)(e?-/cos2a) 

00 00 

+ ^^r^^"^'' (^hal^j Smi{ci, COS a) + ibl^}s'^i{ci,cosa)j sina 

m=l l=m 



(58) 



where, as above, / = 1 for prolate spheroids and / = — 1 for oblate ones. Expressions for 
other factors are given in [13,17]. 

To compare the optical properties of particles of various shape, the cross sections can be 
normalized by the geometric cross section of the equivolume sphere 

C [(ai/6i)2sin2a + cos2«]i/2 



C [{ai/bi)^ cos^ a + sin^ a]^/"^ 



Q for prolate spheroids, (59) 
Q for oblate spheroids. (60) 



vr4 {ai/b,y/^ 
Here, ry is the radius of the sphere with the volume equal to that of a given spheroidal 
particle. This radius can be defined as 



V 



ai bl for prolate spheroids, (61) 
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Vy = a\ hi for oblate spheroids. (62) 

The optical properties of a multilayered confocal spheroid can be found if we put the 
type of spheroid (prolate or oblate), the number of layers A^, complex refractive indices of 
all layers = Uj + kji, the outer aspect ratio ai/hi (ai and hi are the major and minor 
semiaxes), the total particle size parameter, and the relative ratios of volumes of the layers 
fj = Vj/Vtotai- The size parameter may be specified as 

xv = 27rry/A, 

where ry is the radius of a sphere whose volume is equal to that of the spheroid, A the 
wavelength of incident radiation. 

The radial coordinates that define the boundaries of a layered particle, are connected 
with the corresponding semiaxes as 



(l+/)/2 



. 2 

a,- , 



-1/2 

(63) 



The efficiency factors can also be considered as a function of the size parameter 27rai/A given 
by 



The ratio of layer volumes inside the surface j (j = 2, 3, . . . , A^) to the total volume of a 
multilayered particle is determined by using the parameters and C,j 

The aspect ratios of internal layers can be calculated from the volume ratios by using the 
iterative procedure for prolate spheroids 



_ 3 



N 



(0)("-^) + 5^/,[6(^?-l)], (66) 



where n = 1, 2 . . . and the initial value {^j)^^^ = Q-i- For oblate spheroids, the parameter 
can be found by Newton's method 

,in) _ 2[fe-)("-^f + Er^,/z[ei(e? + i)] 

^^'^ ~ 3[(0)("-)]Vl ' ^ ^ 

where n = 1, 2 . . . and the initial value {^jY^^ = 0. 
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6. Numerical Results and Discussion 

6. A. Computational Tests 

The created computer code is based on our codes developed earlier for homogeneous [17] 
and coated spheroids [13]. In calculations of the radial spheroidal functions, we use their 
expansions in terms of the Legendre or Bessel functions, the solution to the corresponding 
differential equation, or Jaffe expansion for prolate functions according to the recommenda- 
tions given in [18,19]. 

The numerical code has been examined by using various tests that include internal con- 
trol (see Table 1), a comparison with the known results for homogeneous and core-mantle 
spheroids and multilayered spheres [20] (Fig. 2) as well as a comparison with the calculations 
for multilayered spheroids based on the EBCM with a spherical basis [1], the SVM with a 
spherical basis [3], and the quasistatic approximation [21]. We also have considered absorbing 
and dielectric particles with a different number of layers and various aspect ratios. 

For non-absorbing particles the efficiency factors for extinction and scattering are known 
to be equal for the same azimuthal index m: Q^^^ = Qi^^ {Q = Q'^^^). Then by increas- 
ing the number of terms A^'^ax in sums for Q^^^ and Qsra, one should obtain a decreasing 



difference between these two factors, i.e., 



/ cxt 
\{m) /^{^) 



; sea 



0, if Nmi,^ — oo. Table 1 shows 



Vext Vs 

the behavior of the efficiency factors in the case of radiation propagating along the rotation 
axis of a spheroid when the sums over m contain only one term, m = 1. A comparison with 
the results presented in [13] demonstrates that the convergence for multilayered spheroids 
resembles that for coated spheroids, i.e, it is not a function of the number of layers A^. The 
convergence is determined by the particle size 27rai/A and is independent of its shape. The 
latter feature makes our solution with a spheroidal basis essentially different from the SVM 
or EBCM approach with spherical basis when convergence quickly degrades for spheroids 
with aspect ratios ai/bi ^ 1.5 — 2 [1,3]. Note that our code allows one to calculate the optical 
properties of spheroids with the size parameters up to xy ~ 15 — 20 including very elongated 
or flattened particles. The time of calculations grows with the increase of layers number as 
t ^ A^^-^, which is much faster as compared to other methods (t ^ A^^.s-s discussion 
in [3]). 

If the particles are nearly spherical, the optical properties of multilayered spheroids and 
spheres should be almost the same. We have considered various absorbing and dielectric 
particles and some results are shown in Fig. 2 where the relative differences in percents 



Qext (sphere) - Cext (spheroid) /vrr^ 
Qcxt (sphere) 



100% (68) 



are given. The values of e are plotted as a function of the size parameter xy for particles 
with the aspect ratio ai/bi = 1.0001 and the volume ratios Vj/Vtotai = 0.33 (j = 1,2,3). 
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The aspect ratios of internal layers are equal to 02/62 ~ 1.00016 and 03/63 ~ 1.00021. The 
wavelike behavior is typical only for dielectric particles. For highly absorbing particles, the 
values of e demonstrate a smooth, monotonous growth with increasing xy- 

6.B. Particles with a Different Number of Layers 

The model of multilayered spheroids gives wide opportunities to investigate both the shape 
and structure effects on the optics of composite particles simultaneously. As one of the first 
applications of the developed model, is our analysis of the idea to represent some composite 
interstellar grains by multilayered particles as suggested in [20]. We consider multilayered 
spheroids with different material layers cyclically changing inside a particle. The particles 
are assumed to be composed of amorphous carbon or silicate with varied volume fraction 
of vacuum. The chosen optical constants for carbon (m = 1.98 + 0.23z) and silicate (m = 
1.68 + 0.03z) correspond to the wavelength A = 0.55 /im. 

The extinction efficiencies of equivolume layered prolate spheroids with ai/61 = 3 are 
compared in Fig. 3. In this case, the aspect ratio of the innermost layer is equal to 03/63 = 
4.91, 09/69 = 8.33 and ai8/6i8 = 11.72 for 3-layered, 9-layered and 18-layered particles, 
respectively. 

The size parameters of compact and porous particles are related as 

.^compact -^compact ( RC\\ 

Xporous = (i_p)l/3 = JyT^^/V^^^^ ^^^^ 

where the particle porosity V if) <V < \) \s introduced as P = Vvac/^totai = 1 — Koiid/Kotai 
and Kac and Koiid are the volume fractions of vacuum and solid material, respectively. 

As in the case of layered spheres (see [20,22]), the scattering characteristics of layered 
spheroids slightly depend on the order of materials and become close to some "average" 
ones, when particles consist of many layers 15 — 20). 

The convergence of the extinction factors seems to be better for oblique incidence and 
oblate particles and larger aspect ratios. Such a behavior is typical for other efficiencies 
(scattering, absorption), albedo and the asymmetry parameter. 

6. C. Particles of Different Porosity 

The difference in the optical properties of compact and porous particles is clearly seen in 
Fig. 4 which shows the size dependence of normalized cross sections as given by Eqs. (59), 
(60). The compact and porous particles have the same mass for the same size parameter. 
This means that variations of the extinction are related to the changes in the particle shape, 
orientation, porosity, and the particle type (prolate or oblate). As follows from Fig. 4, the 
position of the first maximum shifts to larger size parameters with a growth of porosity. For 
very large particles, the normalized cross sections of compact and porous particles cease to 
fluctuate and become rather similar. 
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The role of porosity in dust optics can be properly analyzed by using the normalized cross 
sections 

^(n) C (porous grain) 



C (compact grain of same mass and shape) 
_2/3 Q (porous grain) 



= (l_p)-^/^ -^vi- ' . (70) 

Q (compact grain of same mass and shape) 

The quantity C^'^'^ shows how porosity can increase or decrease the cross section. Such an 
investigation was performed in [22,23] for spherical particles. 

Figures 5, 6 show the extinction cross sections C*^"-* computed for prolate and oblate 
spheroids with porosity V = 0.33 and 0.5. It is seen that the behavior of curves C'-'^^(xcompact) 
is rather complicated. The porosity increases the extinction for spheroids of almost all sizes 
and shapes that are seen pole-on (a = 0°) and decreases the extinction for spheroids that 
are seen edge-on (a = 90°), if xv,compact ^3 — 5. Note that in the last case the curves 
are plotted for the sum of the TM and TE modes for non-polarized incident radiation. For 
very large particles, the normalized cross sections tends to approach to asymptotic values 
C(n) ^ (1 _ p)-2/3 (^ggg Eq. (70)) which are equal to 1.31 and 1.59 if P = 0.33 and 0.5, 
respectively. 

7. Conclusions 

The main results of the paper are as follows: 

1. We have solved the light scattering problem for a confocal multilayered spheroid by 
using the extended boundary condition method (EBCM) with a corresponding spheroidal 
basis. Our recursive solution is based on a special choice of scalar potentials which for any 
next layer can be found by using the potentials of the previous layer, the procedure starting 
from the particle core. 

2. The numerical code has been thoroughly examined by using various tests. They have 
demonstrated that the convergence of the efficiency factors for multilayered spheroids is not 
a function of layers number N and is independent of particle shape. These features make 
our solution with a spheroidal basis essentially different from the SVM or EBCM approach 
with spherical basis. Our code allows one to calculate the optical properties of spheroids with 
the size parameters up to xy = 27rry/A ^15 — 20 {ry is the radius of equivolume sphere) 
including very elongated or flattened particles. 

3. Illustrative calculations show the convergence of the extinction factors to average values 
with a growing number of layers. In the case of the large number of layers, the optical prop- 
erties of layered particles slightly depend on the order of the materials and are determined 
by the volume fraction of the materials. Variations of extinction with a growth of particle 
porosity demonstrate increase of the extinction for spheroids of almost all sizes and shapes 
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that are seen pole-on (a = 0°) and decrease of the extinction for spheroids that are seen 
edge-on (a = 90°), if xy,compact ^3-5. 
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List of Table Captions 

Table 1. Efficiency factors for Extinction Qext and Scattering Qsca for Prolate and Oblate 
Multilayered Spheroids at a = 0°. 
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Table 1. Efficiency factors for Extinction Qcxt and Scattering Qsca for Prolate and Oblate Multilayered Spheroids at a = 0° 







Prolate spheroid 






Oblate 


spheroid 






ai/hi 


= 2 


ai/bi 


= 10 


ai/bi 


= 2 




= 10 


^ 'max 


Qext 


Qsca 


Qcxi 


Qsca 


Qcxt 


Qsca 


Qcx± 


Qsca 


6 


7.3 


7.9 


0.35 


0.41 


2.30 


2.26 


0.26 


0.23 


8 


7.40 


7.36 


0.330 


0.332 


2.409 


2.413 


0.252 


0.251 


10 


7.386 


7.387 


0.3267 


0.3264 


2.4108 


2.4105 


0.2544 


0.2542 


12 


7.38691 


7.38688 


0.32680 


0.32679 


2.410809 


2.410812 


0.25426 


0.25427 


14 


7.3869017 


7.3869022 


0.3268027 


0.3268029 


2.4108093 


2.4108092 


0.254276 


0.254275 


16 


7.38690174 


7.38690174 


0.326802788 


0.326802792 


2.410809320 


2.410809320 


0.25427511 


0.25427512 


18 


7.38690174 


7.38690174 


0.3268027850 


0.3268027850 


2.4108093212 


2.4108093212 


0.2542751277 


0.2542751273 



" the number of layers N = 18, ruj = 1.3 + O.Oi, j = 1,4,7,10,13,16; rrij = 1.5 + O.Oi, j = 2,5,8,11,14,17; mj = 1.7 + O.Oi, j = 3,6,9,12,15,18; 
27rai/A = 5, and V,/Vtotai = 1/18. 



List of Figure Captions 

Fig. 1. Scattering geometry for a prolate spheroid witli tlie confocal layered structure and 
fli/^i = 2. The space is divided into + 1 parts: the outer medium (1), the outermost 
layer (2), ■ ■ -, the core (A^ + 1). The scattered field in the far-field zone is represented in the 
spherical coordinate system {r,-t},ip). 6 is the scattering angle. The origin of the Cartesian 
coordinate system is at the center of the spheroid while the z axis coincides with its axis of 
revolution. The angle of incidence a is the angle between the direction of incidence and the 
z axis in the x - z plane. 

Fig. 2. Percent difference between three-layered spheres and three-layered spheroids e 
defined by Eq. (68): mg = 1.7 + O.Oz, ma = 1.5 + O.Oi, mi = 1.3 + O.Oi, l^j/Kotai = 0.33, 
ai/bi = 1.0001, a = 0°, (•) - prolate spheroids, (o) - oblate spheroids. 

Fig. 3. Size dependence of the extinction efficiency factors for layered prolate spheroids 
with ai/bi = 3. Each particle contains an equal fraction of carbon, silicate, and vacuum (the 
porosity V = 1/3) separated in equivolume confocal layers. The cyclic order of the different 
material layers is indicated (starting from the core). The effect of the increase of the number 
of layers is illustrated. 

Fig. 4. Size dependence of the normalized extinction cross sections for 18-layered prolate 
and oblate spheroids with ai/bi = 3. Particles contain an equal fraction of carbon and silicate 
without vacuum (the porosity V = 0.0) or 50% of vacuum (the porosity V = 0.50). For a 
given value of the size parameter, the compact and porous particles have the same mass. 
The cyclic order of the different material layers is: carbon/vacuum/silicate (starting from 
the core). The effect of the increase of particle porosity and oblique incidence is illustrated. 

Fig. 5. The normalized extinction cross sections (see Eq. (70)) for layered prolate and 
oblate spheroids with ai/bi = 3. For a = 90°, the curves are plotted for the sum of the TM 
and TE modes. The effect of variation of particle type and orientation is illustrated. 

Fig. 6. The normalized extinction cross sections (see Eq. (70)) for layered oblate spheroids. 
The effect of variation of particle shape is illustrated. 
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Incident ray 

Fig. 1. Scattering geometry for a prolate spheroid with the confocal layered 
structure and ai/6i = 2. The space is divided into + 1 parts: the outer 
medium (1), the outermost layer (2), ■ ■ the core (A^+ 1). The scattered field 
in the far-field zone is represented in the spherical coordinate system (r, -d, ip). 
6 is the scattering angle. The origin of the Cartesian coordinate system is at 
the center of the spheroid while the z axis coincides with its axis of revolution. 
The angle of incidence a is the angle between the direction of incidence and 
the z axis in the x - z plane. 
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Fig. 2. Percent difference between three-layered spheres and three- layered 
spheroids e defined by Eq. (68): ma = 1.7+O.Oi, m2 = 1.5+O.Oz, mi = 1.3+O.Oz, 
Vj/Vtota.1 = 0.33, ai/bi = 1.0001, a = 0°, (•) - prolate spheroids, (o) - oblate 
spheroids. 
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Fig. 3. Size dependence of the extinction efficiency factors for layered prolate 
spheroids with ai/6i = 3. Each particle contains an equal fraction of carbon, 
silicate, and vacuum (the porosity V = 1/3) separated in equivolume confocal 
layers. The cyclic order of the different material layers is indicated (starting 
from the core). The effect of the increase of the number of layers is illustrated. 
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Fig. 4. Size dependence of the normalized extinction cross sections for 18- 
layered prolate and oblate spheroids with ai/hi = 3. Particles contain an 
equal fraction of carbon and silicate without vacuum (the porosity V = 0.0) 
or 50% of vacuum (the porosity V = 0.50). For a given value of the size 
parameter, the compact and porous particles have the same mass. The cyclic 
order of the different material layers is: carbon/vacuum/silicate (starting from 
the core). The effect of the increase of particle porosity and oblique incidence 
is illustrated. 
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Fig. 5. The normalized extinction cross sections (see Eq. (70)) for layered 
prolate and oblate spheroids with ai/hi = 3. For a = 90°, the curves are 
plotted for the sum of the TM and TE modes. The effect of variation of 
particle type and orientation is illustrated. 
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Fig. 6. The normalized extinction cross sections (see Eq. (70)) for layered 
oblate spheroids. The effect of variation of particle shape is illustrated. 
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